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I construct initial data for equal-mass irrotational binary black holes using the conformal thin- 
sandwich puncture (CTSP) approach. I locate quasi-circular orbits using the effective-potential 
method, and estimate the location of the innermost stable circular orbit (ISCO). The ISCO predic- 
tion is consistent with results for conformal thin-sandwich data produced using excision techniques. 
These results also show that the ISCOs predicted by the effective-potential and ADM-Komar mass- 
comparison methods agree for conformal thin-sandwich data, just as they did for Bowen-York data. 
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I. INTRODUCTION 

One of the most fervently pursued goals of numerical 
relativity is the simulation of the final orbits, plunge and 
merger of two black holes, and an estimation of the grav- 
itational radiation emitted in the process. The strongest 
gravitational-wave signal is expected to be emitted dur- 
ing the plunge. As such, the most useful numerical sim- 
ulations will describe the binary-black-hole system from 
its last orbit through to the ringdown of the single final 
black hole. It is now possible to evolve binary-black-hole 
spacetimes on an orbital timescale Q, Q , and to extract 
gravitational waveforms from evolutions with even mod- 
est survival times Q- If these waveforms are to be as- 
trophysically realistic, the initial data for the numerical 
simulations must also be realistic. As such, two criteria 
for useful initial data are that (I) they describe two as- 
trophysically realistic black holes in orbit, and (II) this 
should be the last, or almost the last, orbit before the 
black holes plunge together. Due to the circularizing ef- 
fect of gravitational-radiation emission during the black 
holes' inspiral, we expect these orbits to be almost circu- 
lar. 

The initial data used for all of the numerical evolutions 
cited above were those based on the Bowen-York solution 
of the momentum constraint in the conformal transverse- 
traceless (CTT) decomposition of the initial-value equa- 
tions 0, Q ■ These data are suspect on both of the above 
criteria for good initial data. They are known to be as- 
trophysically unrealistic: they are conformally flat, which 
is correct only to second post-Newtonian order even for 
single spinning black holes, and they contain spurious 
gravitational radiation [6| . In addition, predictions of the 
innermost stable circular orbit (ISCO) for these data dis- 
agree with third-order post-Newtonian (3PN) predictions 
by as much as a factor of two 0, and numerical evolu- 
tions of Bowen-York black holes that were predicted to be 
in quasi-circular orbit, are instead found to immediately 
plunge together 0, 0- If we are to obtain meaningful 
results for numerical evolutions, we require more astro- 
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physically realistic initial data, and a better estimate of 
the ISCO. 

A promising alternative to the CTT decomposition is 
York's conformal thin-sandwich (CTS) formulation ja,|9(. 
In all formulations of the initial- value equations, there are 
some pieces of data that we are free to specify as we wish. 
In the CTS formulation, the free data are more closely 
linked to the dynamics of the spacetime than in the CTT 
decomposition, and thus allow us to make physically mo- 
tivated choices. 

The CTS equations have been solved for binary-black- 
hole initial data by a number of groups 10, 11, 12, 13, 
El [H|, and all estimates of the ISCO using CTS data 
are in far better agreement with post-Newtonian results 
[Til fl3L IT^ I than their Bowen-York counterparts. There 
are no published results on the evolution of CTS data, 
but when Bowen-York data are evolved from an initial 
separation greater than the CTS ISCO prediction, they 
are found to orbit [lj, further supporting the CTS ISCO 
over the Bowen- York-data prediction. 

All CTS black-hole solutions to date were found us- 
ing excision techniques: a region around each black hole 
is excised from the computational domain, and inner 
boundary conditions are placed on the resulting exci- 
sion surface. The introduction of excision surfaces is nu- 
merically complex, and one is then presented with the 
problem of choosing appropriate inner boundar y co ndi- 
tions, of which there are several choices [Tol IT3L TTH Il6 | . 
A simpler approach is to use punctures, introduced by 
Brandt and Briigmann E3 to solve the Hamiltonian con- 
straint for Bowen-York data. In the puncture approach, 
divergences at the black-hole locations ( "punctures" ) are 
removed analytically, and the resulting equation(s) can 
be solved numerically over the entire computational do- 
main; no regions need be excised. The puncture ap- 
proach is far easier to implement than excision, and has 
proven useful in constructing Bowen-York binary-black- 
hole initial data Il7l Il8l Il9| , and in numerical evolutions 

eh nam mm. 

In [24[ the puncture approach was extended to the CTS 
equations, and in |25| the conformal thin-sandwich punc- 
ture (CTSP) approach was applied to single boosted- 
black-hole initial-data sets. In this paper I present results 
for &man/-black-hole initial-data sets. In particular, I lo- 
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cate quasi-circular orbits in these data and estimate the 
location of the ISCO. This work provides a useful com- 
parison with excision data results, while avoiding alto- 
gether the issue of excision boundary conditions. It also 
makes it possible for evolution codes that employ punc- 
tures 0,123, to take advantage of CTS data. This 
paper is an extension of the work discussed in |25j | , which 
I will refer to as Paper I. 

There are currently two methods available to lo- 
cate quasi-circular orbits in numerical initial-data sets. 
Cook's effective-potential method [26| identifies orbits of 
particular angular momentum by locating minima in an 
effective potential with respect to the binary's separation, 
in analogy to Newtonian gravity. The mass-comparison 
method of Gourgoulhon, Grandclcmcnt and Bonazzola 
[To| assumes the existence of a helical Killing vector for 
black holes in quasi-circular orbit, and this manifests it- 
self through the equality of the ADM and Komar mass 
estimates. Both the effective-potential | l8U26| and mass- 
comparison [27J methods have been applied to Bowen- 
York data, and are in good agreement in their predictions 
of quasi-circular orbits and the location of the ISCO. For 
CTS data all published work to date has used the mass- 
comparison method to locate orbits 0, 0, 0] . 

It was argued in [24| that it will be difficult (and per- 
haps impossible) to locate orbits in CTSP data using 
the mass-comparison method. However, this argument 
does not imply that orbits do not exist for CTSP data, 
and does not rule out using instead the effective-potential 
method to locate them. As such, in this paper I locate 
quasi-circular orbits in CTSP data using the effective- 
potential method. 

In Section [H] I review the CTSP equations, and the 
procedure for constructing multiple boosted CTSP black 
holes. In Section 1 1 1 1 1 I review the effective-potential 
method. Section IIVI is devoted to numerical methods 
and the setup of the problem for equal-mass irrotational 
binary black holes. Results are presented in Section Ivl 

II. CTSP BLACK-HOLE INITIAL DATA 

The conformal thin-sandwich (CTS) formulation pro- 
ceeds from the 3+1 decomposition of Einstein's equations 
Elll El . In this decomposition spacetime is foliated 
into spacelike slices, and each slice is fully described by 
a spatial metric 7y and an extrinsic curvature -Ky. The 
slices are connected via a shift vector (i % and a lapse func- 
tion N, allowing us to write the spacetime metric as 

ds 2 = -N 2 dt 2 + 7y (da* + (5 l dt){dx 3 + (3 j dt). (1) 

The extrinsic curvature is given by 

Kij = ~ (dt7ij ~ ViPj ~ V, A) , (2) 

where the covariant derivative Vj is with respect to 
the spatial slice. It proves convenient to separate the 



extrinsic curvature into its trace and tracefree parts, 
Kn = Aij + \~tijK. 

In the CTS formulation, the physical quantities 7^, 
Aij, K, N and p % , are all related to conformal back- 
ground quantities via a conformal factor ip: 



lij = i>%j, (3) 

= ip^A.j, (4) 

K = K, (5) 

N = ifj 6 N, (6) 

I? = P l - (?) 



For a good explanation of these conformal weightings, 
and how they relate to the conformal weightings in other 
decompositions, see 0. 

The extra ingredients in the conformal thin-sandwich 
puncture (CTSP) system are the puncture splittings of 
the conformal factor and lapse function. For n black 
holes they are written as 

n 

* = i+E^+ u (8) 

i 
n 

Nif 7 = Nip = 1 + — + v - ( 9 ) 

1 

The quantity r, is the coordinate distance from the ith 
puncture to the field point, rrn parameterizes the mass 
of the ith black hole, and Cj parameterizes the value of 
the lapse at the ith puncture. In the special case of a 
single stationary black hole, the solution u = v = 
and choices mi = M and C\ — ~M provide us with 
the conformal factor and lapse of the Schwarzschild solu- 
tion in isotropic coordinates. In nontrivial examples, the 
motivation for the puncture splittings is that they ana- 
lytically remove the 1/Vj divergences in tp and Nip 7 at 
the locations ("punctures") of the black holes, and allow 
the initial-value equations to be solved numerically over 
the entire computational domain without any need for 
regions around the black holes to be excised. 

In any decomposition of the initial- value equations, 
there are variables that can be freely specified. In the 
CTS formulation the free data are the conformal metric, 
jij, the trace of the extrinsic curvature, K, and their time 
derivatives, Uij = dtjij and dtK. The appearance of «y 
and dtK in the free data is one of the advantages of the 
CTS formulation. Part of the dynamics of the spacetime 
can be freely specified, and this is ideal for constructing 
binary black holes in quasiequilibrium: we simply choose 
Uij = dtK = P.ll3j|. I will make these quasiequilibrium 
choices. For want of a better option, I also choose con- 
formal flatness (7^ = , the flat metric) and maximal 
slicing, K = 0. 

With these choices, the CTSP initial-value equations 
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m vacuum are 



V 2 M 



A L /?-(L/3)^V,lniV 



-l^AvAv, (10) 
0, (11) 



V 2 v 



Nif 7 



(12) 



Equation l|ll)[) is the Hamiltonian constraint, (|llf) is the 
momentum constraint, and (|12|) is the maximal-slicing 
equation. The derivative operators are defined with re- 
spect to the conformal space, i.e., they are flat-space op- 
erators with the choice of conformal flatness. The longi- 
tudinal derivative and Laplacian are defined by 



(L0) ij = V7? J +V J /? l --7 y V fc /? fe , (13) 
A L /f = V^ + iv'V^' (14) 

and conformal flatness gives Rj = 

The conformal extrinsic curvature is constructed from 
the lapse and shift via 



A-- — 
2N 



Mi 



and the choice of ua = reduces this to 
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2N 



ij — ~ {^- d fi)ij ■ 



(15) 



(16) 



When producing nontrivial solutions of the CTSP 
equations, it is necessary for the constants Ci to be pos- 
itive. If at < 0, the conformal lapse will pass through 
zero on some surface(s). A calculation of the conformal 
extrinsic curvature using (|16|) involves division by N, and 
division-by-zero errors will result if N is allowed to pass 
through zero. It was shown in Paper I that different (pos- 
itive) choices of Cj have a negligible effect on the physical 
properties of single boosted-black-hole initial-data sets, 
and the effect of varying Cj on binary-black-hole initial- 
data sets will be considered in Section Ivl 

Note that the CTSP Hamiltonian and momentum con- 
straints can be easily related back to the form of the con- 
straints in the old CTT decomposition [9j . We may solve 
only i|10|) and (fill) , and choose N instead of d^K as a 
piece of free data, removing the need to solve (|12fl . If 
we keep all of the free data choices above, but replace 
dtK = with N = No, a constant, there is an ana- 
lytic solution of the momentum constraint for a boosted 
or spinning black hole: the Bowen-York solution. This 
solution can be used to construct the conformal extrin- 
sic curvature for two boosted black holes, and has been 
used extensively in the past to construct binary-black- 
hole spacetimes. When viewed from the CTS viewpoint, 
we see that Bowen-York data can be considered to satisfy 
one quasiequilibrium requirement, dtjij = 0. The choice 
N = No, however, has no physical motivation. (In- 
deed, it implies an unattractive choice of lapse function, 



N = Nip 6 = NqiJj 6 , which diverges at the punctures.) 
The essential difference between Bowen-York data and 
full CTS (or CTSP) data is due to the addition of d t K 
to the free data, for which a choice can be physically 
motivated. 

For single boosted black holes, it seems that the ad- 
ditional requirement dtK = makes little difference to 
the physical content of the data (see Paper I). However, 
for binary black holes Tichy, Briigmann and Laguna have 
shown |23] that Bowen-York data do not simultaneously 
satisfy both = and dtK = 0. We therefore ex- 
pect CTS and Bowen-York data to differ noticeably in 
the binary-black-hole case. 

Producing boosted-black-hole solutions to the full 
CTSP system (fTU|l . (jTT|) and l(T2*|l also requires a condi- 
tion on the shift vector at each puncture. Boosted black 
holes can be generated by specifying a (non-zero) value 
of the shift vector at each puncture, as outlined in Paper 
I. This shift condition is of the form 



I? = Bl 



(17) 



at the ath puncture, and results in a black hole boosted 
in the B % a j\B % a | direction. The momentum of each black 
hole is parameterized by B % al although that momentum 
is not known analytically and must be calculated numer- 
ically, as was done in Paper I. This procedure allows the 
construction of binary-black-hole initial data for irrota- 
tional black holes — each black hole has non-zero linear 
momentum, but zero angular momentum, with respect 
to a distant observer. It is for equal-mass irrotational 
binary-black-hole spacetimes that I locate quasi-circular 
orbits and estimate the location of the ISCO in Section 
Ivl Before doing that, I will review two methods for lo- 
cating quasi-circular orbits. 



III. QUASI-CIRCULAR ORBITS 

The first method developed to determine quasi-circular 
orbits in binary-black-holc initial-data sets was Cook's 
effective-potential method [2|j. The effective-potential 
method can be motivated by analogy with Newtonian 
gravity [2(| |3(j ■ One can locate circular orbits for a two- 
body system in Newtonian gravity by identifying, for a 
given value of the total angular momentum J of the sys- 
tem, a minimum in the total energy E of the system as 
a function of the separation D of the two bodies. The 
criterion for a circular orbit is then 



dE 
dD 



0. 



(18) 



The angular velocity of the binary system is given by 

dE 



dJ 



(19) 



This procedure leads to Kepler's law for Newtonian or- 
bits. We can use a similar procedure in general relativity 
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if we define an "effective potential" Ef, to take the place of 
the energy E in the Newtonian case, and if we are careful 
to keep the individual masses of the black holes constant 
throughout the process. Cook ,2|| defines an effective 
potential using the total ADM energy of the spacetime, 



E b = E A 



DM 



Mi - Mo. 



(20) 



In the original formulation of the effective-potential 
method, the individual black-hole masses Mi and M2 
were defined in terms of the irreducible mass pfil l3l| . In 
later work Baker [3^ used the puncture ADM masses of 
the black holes. Tichy and Briigmann [2]} found that, 
in the case of Bowen-York puncture data, the puncture- 
ADM and irreducible masses agree to within numeri- 
cal error. In this work I have used the puncture ADM 
masses. 

We also need to choose a suitable measure of the black- 
hole separation, D. This can be the proper distance 
between the apparent horizons of the two black holes 
[Til l26l |. If we are to measure the proper distance, we 
must first find the apparent horizon of each black hole. 
A simpler option is to use, as in this work, the coordinate 
separation between the two punctures [27], |33 • 

Given a binary-black- hole initial-data set, we can iden- 
tify quasi-circular orbits by choosing some total angular 
momentum J, varying the separation of the black holes, 
D (keeping Mi and Mi constant in the process), and 
searching for a minimum in E b , i.e., 



dE b 
3D 



= 0. 



(21) 



This procedure can be repeated for different values of J. 
As J decreases, the black-hole separation of the orbits 
also decreases, but at some point we will find that no 
orbit exists. The last orbit is identified as the innermost 
stable circular orbit (ISCO), which does not exist in New- 
tonian gravity. This method has been applied to both ex- 
cision [2^1 and puncture [Tsl l32l | Bowen-York data, and 
the ISCO predictions from both types of solution are in 
good agreement. 

More recently Grandclement, Gourgoulhon and Bonaz- 
zola OH E3 developed an alternative method to locate 
quasi-circular orbits, which I will refer to as the mass- 
comparison method. The mass-comparison method as- 
sumes the existence of a helical Killing vector in a space- 
time that contains two black holes in quasiequilibrium, 
which is realized through the equality of the ADM and 
Komar estimates of the total energy of the spacetime. 
The mass-comparison method has been applied to ex- 
cision CTS data 0, 0, Q] and puncture Bowen-York 
data HIHI. 

In the case of Bowen-York puncture data, the effective- 
potential and mass-comparison methods agree in their 
predictions of the parameters of quasi-circular orbits and 
the location of the ISCO j2]j. The two methods also 
appear to agree well when applied to CTS excision data 
for corotational equal- mass binaries [34[ • Neither method 
has been applied to CTS-puncture data. 



As was mentioned in the Introduction, it may not 
be possible to apply the mass-comparison method to 
CTSP data. In particular, to locate orbits in the mass- 
comparison method one finds that Cj < [24[, while 
numerical CTSP solutions require Cj > 0, as explained 
in Section [Q] I will instead use the effective-potential 
method. The disadvantage of the effective-potential 
method is that it does not determine a value for the lapse 
parameters Ci in the decomposition ©. If the effective- 
potential method is to give meaningful results, the physi- 
cal parameters of the quasi-circular orbits and ISCO will 
need to be shown to depend only weakly on different 
choices of the parameters Cj. I will discuss this point 
further in Section Ivl 



IV. NUMERICAL SETUP 

I consider equal-mass black-hole binary systems, where 
rrii = = m and ci = Co = c in l(8Jl and 0. The 
CTSP equations {TUJ) - C2> are solved numerically on a 
Cartesian grid, using a multigrid solver that was adapted 
from the BAM_Elliptic solver in the Cactus infrastruc- 
ture [25l |3f|. In Paper I this solver was found to be 
between first- and second-order convergent when solving 
the CTSP equations for boosted black holes. 

The punctures are located on the z-axis at z = ±D/2, 
where D is the coordinate distance between the two punc- 
tures. Each black hole is boosted in the x-direction by 
specifying 



(22) 



at its respective puncture. In the single-black-hole results 
presented in Paper I, the puncture was located at the 
origin; condition 1)22(1 could be imposed on every level 
of the multigrid hierarchy because the origin existed on 
every level. In the binary-black-hole case the puncture is 
placed on a point that exists on the finest grid, but may 
not exist on any of the coarser subgrids. Condition i|22|) 
is applied on the finest grid, and on each of the coarser 
grids the value of the shift vector is specified at one of 
the points nearest to the puncture, and given a value 
such that first-order interpolation would give the correct 
value (T-B) at the puncture. 

The CTSP equations are solved on only one octant of 
the Cartesian grid (i.e., for only positive x, y and 2), by 
using the symmetries of the CTSP variables consistent 
with a binary system with non-zero angular momentum, 
shown in Table [J These symmetries are found by not- 
ing that as r — -> 00, the CTS momentum constraint 
approaches the form of the CTT momentum constraint, 



Al/3* = 0, 



(23) 



for which we have the Bowen-York solutions for a boosted 
or spinning black hole. Far from the binary, I assume that 
the solution of the momentum constraint will approach 
that for a single spinning black hole |36|. and that the 
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x = 


v = o 


2 = 


r — > oo 


shift vector are 


u 


even 


even 


even 


1/r 




V 


even 


even 


even 


1/r 


N l di{r 2 (3 l ) 




even 


even 


odd 


z/r 3 


py 


odd 


odd 


odd 


1/r 3 


for /? x and f3 z , and 


P z 


odd 


even 


even 


x/r 3 



TABLE I: Coordinate plane symmetries and outer boundary 
conditions for CTSP variables for an equal-mass irrotational 
binary-black-hole system with angular momentum in the y- 
direction. 



shift vector has the angular and radial dependence of the 
vector potential for a single spinning Bowen-York black 
hole, 



ijk 



rijJk 



(24) 



where J % is the direction of the total angular momentum 
of the system and n l is the radial normal vector directed 
away from the origin. The same symmetries were used 
by Baumgarte, et. al. [37j to solve the conformal thin- 
sandwich equations for binary neutron stars. For the 
numerical setup I have described, ,P — (0, 1,0). 

In the multigrid solver, careful treatment of the sym- 
metry planes was required whenever a component of the 
shift vector was odd across that plane. The odd symme- 
try requirement is that the function has equal and oppo- 
site values on either side of the symmetry plane, and be 
zero on the plane. However, at intermediate steps in the 
solution procedure the function is unlikely to be exactly 
zero at the points on the symmetry plane. I have found 
that an improvement in efficiency of the code by many or- 
ders of magnitude can be achieved by setting the function 
to zero on the symmetry plane before any prolongation 
operation (which interpolates the function from a coarse 
to a fine grid), but not at any other time. 

At the outer boundary, the functions u and v fall off 
as 1/r. The behavior of the shift vector at the outer 
boundary is required to be consistent with i|24|) . For the 
setup described, the condition l|24[l applies to only two 
components of the shift vector, (3 X and f3 z ; the 1/r 2 con- 
tribution to f3 v is zero. As such, the code requires that 
(3 V fall off as 1 /r 3 . The outer boundary behavior of the 
CTSP variables is shown in Table [I] 

To enforce the asymptotic behavior shown in Table [I] 
the code applies Robin outer-boundary conditions to all 
of the CTSP variables. The conditions applied to u and 
v are 



rfdiiru) = 0, 
N l d l (rv) = 0, 



(25) 
(26) 



where N l is the unit vector perpendicular to the bound- 
ary. The conditions applied to the components of the 



N^iir 3 ^) = 



(27) 



(28) 



for py. 

For a given solution of the CTSP equations with the 
shift condition (12211 and outer boundary conditions l|25[) . 
(|2(i|l . (|27l) and 128|) . we can calculate a number of physical 
quantities. The total angular momentum of the space- 
time is given by 0, |2(J 



Ji 



£ »jfc 
8tt 



x J A M d 2 S t 



(29) 



The surface integral is computed at the outer boundary 
of the computational grid. The integral was constructed 
using global Killing vectors of the conformal space (which 
are asymptotic Killing vectors of the physical space) and 
can in fact be computed at any radius, so long as the 
surface surrounds the puncture. For the setup we have 
described, the symmetries will enforce J x = J z = 0. 

The total ADM mass of the spacetime is given by 0, 
381 



n 

M ADM ^m, + -j iP^AijA^dV. 



(30) 



The volume integral in l|30|) is evaluated over the entire 
numerical grid. The contribution from outside the nu- 
merical grid is estimated with the aid of the outer bound- 
ary conditions, using a technique described in Paper I. 

There is no unique measure of the mass of each indi- 
vidual black hole. I have chosen to use the ADM mass 
computed on the second hypersurface of each black hole, 
sometimes called the bare mass 0, |^| • I will refer to 
it as the puncture ADM mass, because it can be easily 
computed at each black hole's puncture. The puncture 
ADM mass, M i; is 113 



Mi = mi ( 1 + u\ 



m 2 \ 
2DJ 



(31) 



where u\ is the value of the function u at the first punc- 
ture, and M2 is found by interchanging the subscripts "1" 
and "2" . Note that due to the symmetries of the equal- 
mass binary case, M\ = Mi. Define the total black- 
hole mass as M — M\ + M2 and the reduced mass as 
fx = M x M 2 /{Mi + M 2 ), which are M = 2 and fi = 1/2 
in the results presented below. 

When producing lines of constant J and M in the 
effective-potential method, the constants B in l|22(l and m 
in (JSJ) are varied until J and M are equal to the required 
values, with a combined error of less than 3 x 10~ 6 . 
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D 



J 



Eadz 



Method 



J 



0.25 1.625 3.1175 -0.0774 1.961326 0.1144 

0.5 1.625 3.1175 -0.0773 1.961342 0.1137 

1.0 1.625 3.1150 -0.0775 1.961228 0.1118 

2.0 1.625 3.1125 -0.0777 1.961169 0.1133 

TABLE II: ISCO results for various values of the lapse pa- 
rameter c, with resolution h — 0.125m and outer boundary 
at 32m. 



V. RESULTS 

In the first set of results presented here, the lapse pa- 
rameter is c = 1.0. Quasi-circular orbits were found for 
numerical solutions with a resolution of h = 0.0625M and 
an outer boundary at 16M. The innermost stable circu- 
lar orbit was identified at J = J/(/iM) = 3.115 ± 0.005, 
D = D/2M = 1.625 ± 0.125, and E b = E b /n = 
-0.0775 ±0.002. The orbital angular speed at the ISCO 
was O = MSI = 0.11 ± 0.01. 

The uncertainties in these quantities were found by 
locating the ISCO for solutions with the same resolu- 
tion but a more distant outer boundary (h = 0.0625M, 
outer boundary at 2AM). There is also uncertainty due 
to the freedom in the lapse parameter c. As discussed 
in Section IlIII the lapse parameter c is undetermined by 
the effective-potential method. I investigated the depen- 
dence of the location of the ISCO on this parameter for 
the values c = 0.25,0.5,1.0, and 2.0. Note that the el- 
liptic CTSP system becomes more difficult to solve as c 
decreases, and I could not find solutions for c = 0. This 
is unfortunate, because when c = the lapse would be 
zero at the black-hole punctures, and constitute a fully 
"pre-collapsed lapse", a potentially desirable choice for 
numerical evolutions 0, HO, |4£J |4l| . 

The ISCO parameters for different choices of c are 
shown in Table [H] It is clear that the value of c has 
only a small effect on the ISCO location, and the values 
remain well within the uncertainties quoted above. 

Table Hill compares ISCO results for irrotational bina- 
ries, obtained by different methods. The "BY-Puncture" 
results at the end of the table are those obtained by 
Baumgarte for Bowen-York puncture data. These 
results are not close to any of the post-Newtonian re- 
sults, from either effective one- body (EOB) Q or stan- 
dard post-Newtonian methods. The CTSP results 
in this paper and the CTS excision results of Cook and 
Pfeiffer hj agree in their estimate of fl at the ISCO. 
They do not agree on J and E b , but are very close. Both 
sets of results are much closer to the post-Newtonian es- 
timates than the Bowen- York-puncture results — in fact, 
they are between the 2PN and 3PN results from the 
"standard" and effective one-body methods. The vari- 
ous ISCO values of E b and f2 from different methods are 
also shown in Figure ^ 



CTSP 3.115 
CTS-Excision 3.06 



0.0775 0.112 
■0.0724 0.101 



2PN (standard) 3.116 
3PN (standard) 3.144 



0.0796 0.137 
0.0772 0.129 



2PN (EOB) 3.408 
3PN (EOB) 3.28 



0.0576 0.0732 
■0.0668 0.0882 



BY-Puncture 2.95 -0.092 0.18 

TABLE III: ISCO results for irrotational equal-mass binaries 
from different methods. The CTS-Excision results are typical 
values from Cook and Pfeiffer [l3l ]. The "PN (standard)" re- 
sults are from Blanchet |42| using a standard post-Newtonian 
expansion. The "EOB" results are from effective one-body 
post-Newtonian methods 0. The BY-puncture results are 
for Bowen-York puncture data [T^|. 



-0.070 



: 2PN (standard) 
■ 3PN (standard) 

2PN (EOB) 
I3PN (EOB) 
f BY-puncture 
I CTS-Excison 
sCTSP 

0.1 



FIG. 1: ISCO results for equal-mass irrotational binaries. The 
methods are as described in Table [m] The error bars for the 
CTSP result are based on the uncertainties quoted in the text, 
with c = 1.0. 



VI. DISCUSSION 

I have numerically solved the CTSP equations to con- 
struct initial-data sets for irrotational binary-black-holc 
spacetimes. Contrary to the implications of [24(. it is 
possible to locate quasi-circular orbits in CTSP data, by 
using the effective-potential method instead of the mass- 
comparison method. Identification of quasi-circular or- 
bits allows an estimate of the innermost stable circular 
orbit (ISCO). The physical parameters of the ISCO are 
close to those found by Cook and Pfeiffer 0] > wn0 solved 
the CTS equations with excision techniques and used the 
mass-comparison method to locate quasi-circular orbits. 
This suggests a great deal of robustness of the results: 
essentially the same CTS ISCO was found using two dif- 
ferent solution procedures (punctures and excision), and 
two different methods of determining quasi-circular or- 
bits (the effective-potential and mass-comparison meth- 
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ods). These results also suggest that, as in the case of 
Bowen-York data j2?|], the effective-potential and mass- 
comparison methods agree in their determination of 
quasi-circular orbits. 

Although the mass-comparison method was not used 
to locate quasi-circular orbits in the CTSP data, it is 
natural to wonder how close the ADM and Komar mass 
estimates are for these data. Recall that in the mass- 
comparison method equality of the two mass estimates is 
considered to denote a quasi-equilibrium spacetime |10| . 
One immediately finds that in the CTSP approach, the 
ADM and Komar masses cannot be equal. This can be 
seen by writing out the Komar mass as was done for the 
ADM mass in . The ADM mass for an equal-mass 
binary can be written 

M ADM = 2m + I 1 , (32) 

where I\ is the integral in l|30[) . and m = m\ = m 2 is 
the mass parameter. The Komar mass can be written 

M K = m-c + h/2 + I 2 , (33) 

where I2 is similar to I\. In particular, I\ and 1% are both 
positive, and small compared to the mass parameter m. 
If Madm and Mk are to agree, we need roughly c = — m. 
Using c > 0, as in the CTSP procedure, Madm and Mk 
will always disagree. They will be closest for the smallest 
choice of c (which is 0.25 in the results reported here), 
but still differ by more than a factor of two. 



The results suggest that the physical content of the 
data changes little as c is changed, while the Komar mass 
will change a great deal. One could argue that the data 
with the closest ADM and Komar masses (i.e., when c is 
smallest) represent the data closest to quasiequilibrium. 
However, it is not clear how one could quantify what 
"close" means in this case. 

It is now possible to construct CTS binary-black-holc 
initial data for use in codes that use punctures. In punc- 
ture evolutions the lapse quickly collapses to zero at the 
punctures |20l .23], and codes can become unstable if 
the lapse is below zero anywhere on the numerical grid. 
CTSP data are well suited to both of these requirements. 
For suitable choices of the parameters Cj, the lapse is par- 
tially "pre-collapsed" at the punctures, and the lapse by 
construction must always be positive. Evolution of these 
data sets will be the subject of further research. 
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